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ABSTRACT 

We show that a survey of equations of state and observations of X-ray pulsations from SAX J1808.4- 
3658 give 2.27 Mq as the upper Hmit of the compact star mass. The corresponding upper hmit of the 
radius comes out to be 9.73 km. We also do a probabilistic study to estimate the lower limit of the 
mass of the compact star. Such a limit puts useful constraints on equations of state. We also discuss 
the implications of the upper mass limit for the the evolutionary history of the source, as well as the 
detection of it in radio frequencies. We envisage that the possible observation of radio-eclipse may be 
able to rule out several soft equation of state models, by setting a moderately high value for the lower 
limit of inclination angle. 

Subject headings: accretion, accretion disks — binaries: close — equation of state — pulsars: 
individual (SAX J1808.4-3658) — X-rays: stars 
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1. INTRODUCTION 

The discovery of millisecond X-ray pulsations (period 
T = 2.49 ms; Wijnands & van der Khs 1998) in the tran- 
sient X-ray burster SAX J1808. 4-3658 confirmed the spec- 
ulation that LMXBs are progenitors of millisecond pul- 
sars (Bhattacharya & van den Heuvel 1991). The orbital 
period (Porb = 2.01 hr) and the pulsar mass function 
(/i = 3.7789 X 10~^) of this source were observationally 
determined by Chakrabarty & Morgan (1998). These give 
valuable information about the masses (of both the pri- 
mary and the secondary) and the inclination angle. For 
example, the value of Porb uniquely determines the mass of 
a Roche lobe filling low mass star with known mass-radius 
relation. 

It has been recently proposed that the compact star in 
SAX J1808. 4-3658 is a strange star (SS) and not a neu- 
tron star (NS) (Li et al. 1999). Such a speculation, if 
confirmed, will prove that the so called strange matter hy- 
pothesis (Witten 1984) is correct. According to this hy- 
pothesis, strange quark matter (made entirely of decon- 
fined u, d and s quark) could be the true ground state of 
strongly interacting matter rather than ^^Fe. This is an 
important problem of the fundamental physics. To resolve 
it, we need to constrain the equations of state (EOS) for 
this compact star very effectively. 

In this Letter, we estimate the upper limits of the mass 
and the radius of the compact star in SAX J1808. 4-3658. 
We also discuss the possible ways to estimate lower mass 
limit. 



2. UPPER LIMITS TO MASS AND RADIUS 

We estimate the upper limits of the mass and the ra- 
dius of the compact star in SAX J1808. 4-3658 by using 
the basic requirements for X-ray pulsations. Here, upto 
eqn. 5, we follow the same method as described in Li et 
al. (1999). To explain it, we first define the corotation 
radius (Pco) and the magnetospheric radius (Pmag)- They 



are given by (see Burderi & King 1998) 
Pco = 1-5 X lOSn^^T^/^ 



and 
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(1) 



(2) 



Pmag = 1.9 X 10>Ai26 

where mi is the compact star mass in units of solar mass, 
T is the compact star spin period in units of milliseconds, 
(f) is the ratio between the magnetospheric radius and the 
Alfven radius, /i26 is the compact star magnetic moment 
in units of 10^^ G cm^ and Mn is the accretion rate in 
units of 10^^ g s^^. Pco and Pmag are given in cm. In 
this Letter, we assume that (j) is almost independent of the 
accretion rate (Burderi & King 1998). 

The requirements for X-ray pulsations (if there is no 
"intrinsic" pulse mechanism) and the presence of accre- 
tion flow (that is not centrifugally inhibited) give (see Li 
et al. 1999) 



Pi < Pmag(Afniax) < Pmag (^^^min) < Pc 



(3) 



where M^-m and Mmax give the range of the accretion rate 
in which X-ray pulsations in SAX J1808. 4-3658 were ob- 
served. From eqns. 1-3 we get (Li et al. 1999) 



Pi < 27.1 



T 



(4) 



P„,i„' ^ 2.49ms' 
where Pmax and Pmin are maximum and minimum values 
of measured X-ray fluxes respectively. It is to be noted 
that here the pulsar magnetic fleld is assumed to be dipo- 
lar. In writing eqn. 4, we also assume that M is pro- 
portional to the observed flux F for all accretion rates. 
This is justified by the fact that the X-ray spectrum of 
SAX J1808. 4-3658 was remarkably stable (Gilfanov et al. 
1998), when the X-ray luminosity varied by a factor of 
~ 100 during the 1998 April/May outburst. During this 
period, in 2—30 keV band, the maximum observed flux was 
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around 3 x 10~^ erg cm~^ s~^, while the flux value dropped 
to around 2 x 10~^^ erg cm^^ when the pulse signal be- 
came barely detectable (Cui, Morgan, & Titarchuk 1998). 
Adopting maximum to minimum flux ratio as 100 (after 
Li et al. 1999), we get from eqn. 4 

Ri < TAOmY^ km (5) 
Eqn. 5 gives the maximum value of Ri, if the maximum 
value of mi is known. To calculate mi, max we flrst rewrite 
the eqn. 5 in the following form. 

mi < 11. 19x1;^/^ (6) 
where xi is the dimensionless radius to mass ratio of the 
compact star. We can compute mi, max from eqn. 6, if 
the minimum value of Xi is known. To choose the value 
of a;i,min, we survey about 20 EOS (that include both SS 
and NS) and examine the value of xi corresponding to 
the maximum possible mass for a given EOS. For both 
SS and NS, we choose EOS of widely varying stiffness pa- 
rameters, which guarantees our results to be of suflicient 
generality. This is reflected by the wide range of maxi- 
mum possible mass values given in Table 1, where we have 
listed 13 representative EOS. From Table 1 and Figure 1, 
we notice that the xi values for all the EOS are confined 
to the range 2.98 — 4.34, with 11 (out of 13) points clus- 
tering in 3.3 — 3.7. To illustrate this, we draw a vertical 
line in Figure 1, corresponding to xi = 2.9. As none of 
the EOS points falls to the left of this line, we take 2.9 as 
the lower limit of xi. Such a conclusion is very general, as 
it is valid for the whole range of existing EOS. This gives 
2.27 (the crossing point of the vertical line and the curve 
in Figure 1) as the upper limit of mi from eqn. 6. The 
corresponding upper limit of Ri comes out to be 9.73 km 
from eqn. 5. 

It is to be noted that for some SS EOS, xi value may be 
less than 2.9 for lower values of masses (i.e., less than the 
maximum possible mass). But as wc use the lower limit 
of xi to estimate the maximum possible value of mi, it 
is justified to take 2.9 as the minimum possible value of 
xi. An EOS model (that may be put forward in future) 
with xi (corresponding to the maximum possible mass) 
less than 2.9, will give a higher value of mi, max than 2.27. 
However, such an unusual EOS is highly improbable. We 
also point out that if we take into account the rotation of 
the compact star, the lower limit of xi will increase, re- 
sulting in a decrease of mi, max- Therefore we can say that 
2.27 may be the firm upper limit of mi. 

For the sake of completeness and to give more credibility 
to our work, we calculate mi, max with less constraining val- 
ues of xi. For this purpose, we take a;i,,jii,i ~ 2.25, which 
is the absolute lower limit (for compact star) of xi (Wein- 
berg, S. 1972). This hmit, which is independent of EOS 
and depends only on the structure of the relativistic equa- 
tions for hydrostatic equilibrium, gives mi, max = 3.32 and 
^i,max — 11.04 km. Therefore 3.32 is the absolute upper 
limit of mi. Another value of xi,min(= 2.56) was derived 
by Bondi, H. (1964), under the reasonable assumptions 
concerning the EOS, i.e., e > 0, p > and dp/de < 1 
(where p is pressure and e is energy density). This value 
of a;i,min implies mi, max = 2.73 and i?i,max = 10.34 km. 
Therefore we see that the value of i?i,max is not very sen- 
sitive to the chosen value of a;i,min- 



3. LOWER MASS LIMIT 

Here we estimate the probability (Pmin) of a possible 
compact star mass (m,i) to be the lower limit of mass. Wc 
do it using the procedure " Random distribution of orbital 
inclinations" for measuring neutron star mass, mentioned 
in Thorsett & Chakrabarty (1999). Because of the ab- 
sence of sufficient observational data, here wc can not fol- 
low any well-established statistical method. For example, 
the measured value of a single post-Keplerian (PK) pa- 
rameter (Taylor, J.H. 1992) (with additional assumptions, 
such as a uniform prior likelihood for orbital orientations 
with respect to the observer) can be used to make strong 
statements about the posterior distribution of the masses 
(Thorsett & Chakrabarty 1999). But none of these param- 
eters could be measured for SAX J1808.4-3658. Therefore, 
our results basically depend on the a priori probability of 
observing the source with a given inclination angle (i). 

To explain the method, we first rewrite the well known 
expression for the pulsar mass function (/i) in the follow- 
ing way. 

smi = // ^ — (7) 

m2 

where m2 is the mass of the companion star in units of 
solar mass. For a main sequence companion that fills its 
Roche lobe, m2 = 0.22 (corresponding to Porb = 2.01 hr). 
As a result, the lower limit of i (i.e., Win) comes out to 
be 3° from eqn. 7 (using mi,min = 0; the absolute lower 
limit). However Chakrabarty & Morgan (1998) have ar- 
gued that m2<0.1 (because the companion is bloated by 
irradiation). Therefore, we take 0.1 as the upper limit 
of m2, that corresponds to Zmin = 4". The absence of a 
deep eclipse indicates that for a Roche lobe filling com- 
panion, imax = 82° (Chakrabarty & Morgan 1998). We 
set m2,min = 0.05, which is a possible companion mass ac- 
cording to Chakrabarty & Morgan (1998). We also take 
two other values of m2,min for the purpose of illustration. 

With all these limiting values, we calculate Pmin in the 
following way. Given a value of mi, we compute the al- 
lowed range of i (i.e., ia,min and ia,max) from eqn. 7, for the 
chosen range of m2 (i.e., m2,min < m2 < m2,max)- How- 
ever, for ia,max > «max, WC take ia,max = «max- Similarly, 
we do not consider any value of i less than imin- Now we 
argue (after Chakrabarty & Morgan 1998) that in statisti- 
cal calculations it is useful to assume that binary orbits are 
randomly oriented with respect to the line of sight (sec also 
Thorsett & Chakrabarty 1999). The differential distribu- 
tion of inclinations is then proportional to sini. This gives 
the a priori probability of observing a system with i in the 

range ^a,min — ^ ^ ^a.max ^S P — (cOS Za,min COS Za, max)- 

Therefore P should be the probability of the chosen mi 
for being the actual compact star mass. 

We calculate P for many mi values (at regular inter- 
val) in the range mi,min < mi < mi, max- Then Pmin is 
calculated by the formula 

V" P 

^min p V°/ 

where P is calculated at n number of mi points and j is the 
index number of the mi value at which Pmin is required. 
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In figure 2, we plot Pmin against mi for tliree m2,min 
values (0.04, 0.05 and 0.08). For TO2,min = 0.05, we see 
that the minimum value of mi is 0.70 with 95% probabil- 
ity. Therefore, although we start with zero as the lower 
limit of mi, we get a large value for mi with a very 
high probability. This shows that the probabilistic method 
is very effective in estimating the value of mi^min- If we 
take mimax = 3.32 (i.e., the absolute limit), the minimum 
value of mi comes out to be 0.90 with 95% probability, 
which shows that this method is sensitive to the assumed 
value of mi, max (a less constrained value of the upper 
limit of mi gives a higher value of mi^min with the same 
probability). A considerable increase in either m2.max or 
«max does not change our result much. For example, if 
?Ti2,max = 0.22, the 95% probabilistic minimum value of 
mi is 0.68. However our result changes with the change of 

m2,min for m2,min < 0.04. 

We also point out that for a given value of m2.min, if 
*min is greater than a certain value, it can be seen from 
eqn. 7 that every ij^in will correspond to a minimum 
possible value of mi. Therefore if we can observation- 
ally constrain i from the lower side, the value of mi.min 
can be predicted more accurately (as it will not depend on 
the probabilistic study). For example, the detailed mod- 
eling of the optical companion's multiband photometry 
during outburst with a simple X-ray heated disk model 
suggests that cosz < 0.45 (Wang et al. 2001, in prepa- 
ration; Bildstcn & Chakrabarty 2001) for SAX J1808.4- 
3658. This implies i > 63° and hence mi^min = 1-48 (using 
m2,min = 0.05). 

4. DISCUSSIONS 

In this Letter, we have estimated the upper limits of the 
mass and the radius of the compact star in SAX J1808.4- 
3658. Li et al. (1999) have concluded that a narrow region 
in TOi — Ri space will be allowed for this star. The upper 
boundary of the mass will constrain this region effectively. 
It can also give the upper limit of i (from eqn. 7), if m2,min 
is known by an independent measurement. Alternatively, 
"T'l.max gives the upper limit of m2, for a known value of 
^min- For example, in^n = 63° gives m2,max = 0.066(0.085) 
for mi,max = 2.27(3.32). 

As we have mentioned in section 2, in this work, the pul- 
sar magnetic field is assumed to be primarily dipolar. If 
the field has more complicated structure, the Rmag — M re- 
lation will be changed, resulting in the modification of eqn. 
4. This will lead to the change in eqn. 6, and hence our 
calculated value of wi^max (and i?i,max) will be modified. 
However, Li et al. (1999) have argued that the accretion 
flow around the compact star is dominated by a central 
dipole field, which gives credibility to our results. 

Corresponding to every EOS, there exists a maximum 
possible mass. Therefore a lower limit of mi is very im- 
portant in constraining EOS and hence in understanding 
the properties of matter at very high density compact star 
core. The possibility of this candidate to be a strange star 



can also be checked more effectively. Using our figure 2, 
we can predict the probability with which a certain value 
of mi will be the lower limit. For example, mi = 1.41 
(the maximum possible mass for our model Y) will be the 
lowest possible mass with 72% probability (from curve 1 
of Figure 2). However, it is to be kept in mind that such 
a probabilistic study may be valid, if binary inclination 
angles are distributed randomly. 

If the orbital evolution of SAX J1808.4-3658 is driven by 
only gravitational wave radiation, then the rate of change 
of the orbital period will be given by (Ergma & Antipova 
1999) 

Porb = -1.72 X 10-^(27r/Porb)^/^m2mi(m,i + m2)-^'^) 
This implies 2.30 x 10~^^ as the maximum possible (ab- 
solute) value of Porb (for mi, max = 2.27, m2,max = 0.1 and 
Porb = 7249 s). Chakrabarty & Morgan (1998) have sug- 
gested that Porb can be measured, if the source remains in 
the X-ray bright state for long enough (or if the pulsations 
remain detectable in quiescence). If in future, such a mea- 
surement yields the value of the orbital period decay rate 
greater than 2.30 x lO'^^ (2.99 x 10"" for mi,max = 3.32), 
then we can conclude with certain confidence that the or- 
bital evolution of SAX J1808. 4-3658 is significantly driven 
by magnetic braking. This will give support to the evolu- 
tionary model of Ergma & Antipova (1999) and in general 
will be very important for learning about the prehistory 
of the system. A better understanding of the criterion for 
magnetic braking will also be possible. 

It has been proposed that SAX J1808. 4-3658 may 
emerge as a radio pulsar during the X- ray quiescence 
(Chakrabarty & Morgan 1998). Ergma & Antipova (1999) 
have calculated that for A < 3 cm, it may be possible 
to observe radio emission from this source. However our 
limits of mass values give a slightly higher (3.8 cm for 
m\ max = 2.27 and 4.5 cm for mi max = 3.32) upper limit 
for A. 

As we have already mentioned in section 3, a moder- 
ately high value of imin will give a lower limit of mi (with- 
out any probabilistic study). This will be very impor- 
tant for constraining EOS more decisively. For example, if 
imin = 63° (corresponds to mi^min = 1.48, given in section 
3), our EOS models SSI, SS2 and Y will be unfavoured (see 
Tii.max-column of Table 1). According to Chakrabarty & 
Morgan (1998), a deeper eclipse might be observed for the 
less penetrating radio emission, providing a strong con- 
straint on the value of i. Therefore, we expect that, the 
value of Zmin (determined by this method) may be able to 
rule out several soft EOS models in future. 

We acknowledge Dipankar Bhattacharya for reading the 
manuscript and giving valuable suggestions. We thank Xi- 
angdong Li for sending useful comments and Pijush Bhat- 
tacharjce for encouragement. We also thank the referee 
for giving suggestions to improve the quality of the paper. 
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Fig. 1. — nil vs. x'l plot (see the text) for compact star. The solid curve indicates the upper bound of the mass according to eqn. 6. The 
asterisks are for different EOS models listed in Table 1. The vertical line corresponds to xi = 2.9. 



6 S. Bhattacharyya 

Table 1 

List of 13 EOS (both SS and NS) of widely varying stiffness parameters and their references. The values of 

RELEVANT PROPERTIES (SEE THE TEXt) ARE ALSO GIVEN. 



EOS label 


compact star 


Reference 


^l,max 


Ri (km) 




SS2 


SS 


Dey et al. (1998) 


1.32 


6.53 


3.34 


SSI 


SS 


Dey et al. (1998) 


1.44 


7.07 


3.32 


Bgo 


SS 


Farhi & Jaffe (1984), S = 90 MeV/fm^ m, = 


1.60 


8.74 


3.69 


Beo 


SS 


Farhi & Jaffe (1984), B = 60 MeV/W, m, = 


1.96 


10.71 


3.70 


Y 


NS 


Pandliaripande (1971b), liyperonic matter 


1.41 


7.10 


3.39 


B 


NS 


Baldo, Bombaci, & Burgio (1997), nuclear matter 


1.79 


9.64 


3.64 


W 


NS 


Walecka (1974), neutron matter 


2.28 


11.22 


3.32 


SBD 


NS 


Sahu, Basu, & Datta (1993), nuclear matter 


2.59 


14.08 


3.68 


A 


NS 


Pandharipande (1971a), Reid soft core 


1.66 


8.37 


3.42 


AU 


NS 


Wiringa, Fiks, & Fabrocini (1988), AV14 + UVII 


2.13 


9.41 


2.98 


FPS 


NS 


Lorenz, Ravenliall, & Petliick (1993), UV14 + TNI 


1.80 


9.28 


3.48 


L 


NS 


Pandharipande & Smith (1975b), mean field 


2.70 


13.70 


3.43 


M 


NS 


Pandharipande & Smith (1975a), tensor interaction 


1.81 


11.60 


4.34 
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Fig. 2. — Pmin vs. mi plot with the following parameter values : mi^min = 0, mi, max = 2.27, m2,max = 0.1, imin = 4° and imax = 82°. 
Curve 1 is for m2,min = 0.05, curve 2 for m2,min = 0.08 and curve 3 for m2,min = 0.04. 



